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SUMMARY 


Two  aspects  of  power  flow  associated  with  electromagnetic  waves  in  plane- 
stratified,  dispersive,  anisotropic  media  that  are  also  lossless  and  linear  are  considered. 
One  aspect  is  the  relation  between  group  velocity  and  the  velocity  of  energy  transport  of 
surface  waves  in  such  media.  It  is  shown  that  the  group  velocity  of  surface  waves  is 
equal  to  the  ratio  of  the  real  part  of  the  complex  Poynting  vector,  integrated  over  the 
coordinate  of  stratification,  to  the  corresponding  integral  of  the  stored  energy  density. 

Tha  second  aspect  is  the  relation  between  the  dyadic  surface  impedance  representing 
either  a  slab  cf  plane-stratified  medium  above  a  perfectly  conducting  plane  or  a  semi- 
infinite  region,  the  latter  for  the  case  of  evanescent  fields,  and  the  pewer  flow  in  the 
respective  structures.  The  significance  of  the  surface  impedance  and  power  relations 
for  surface  waves  is  discussed. 
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INTRODUCTION 


In  this  paper  two  aspects  of  the  power  flow  associated  with  electromagnetic 
waves  in  plane-stratified,  anisotropic,  dispersive  media  and  their  application  to  surface 
wave  propagation  are  considered.  The  first  aspect  is  that  of  the  relation  between  the 
gro^p  velocity  and  the  velocity  of  energy  transport  of  surface  waves;  the  second  is  the 
relation  between  a  dyadic  surface  impedance  and  the  power  flow  and  stored  energy  in  the 
structure  it  represents. 

It  is  well  known  that  monochromatic  plane  electromagnetic  waves  in  a  homogeneous, 

dispersive,  anisotropic  medium  that  is  also  lossless  and  linear,  e.  g.  ,  the  ionosphere  for 

small- signal  propagation,  carry  power  in  the  direction  of  the  normal  to  the  plane  wave 

dispersion  surface.  Specifically-  the  velocity  of  energy  transport  of  plane  waves  in  such 

a  medium  is  equal  to  the  group  velocity,  that  is,  the  gradient  in  the  wave  number  space 

of  the  frequency.  ^  This  relation  between  the  group  velocity  and  the  velocity  of  energy 

transport  finds  an  important  application  in  the  ray  interpretation  of  the  far  fields  radiated 

(3) 

by  sources  in  the  presence  of  homogeneous  anisotropic  media. 

Ic  is  shown  here  that  an  analogous  relation  involving  the  group  velocity  holds  for 
the  case  of  surface  waves  in  plane- stratified,  dispersive,  anisotropic  media  that  are  also 
lossless  and  linear,  in  that  the  group  velocity  of  surface  waver  that  can  propagate  in  such 
a  medium  may  be  interpreted  as  the  surface  wave  energy  transport  velocity.  For  such 
surface  waves,  the  direction  as  well  as  the  magnitude  of  the  real  part  £  of  the  complex 
Poyntxng  vector  is,  in  general,  a  function  of  z,  the  coordinate  in  the  direction  of  strat¬ 
ification.  (An  example  of  this  dependence  is  described  in  Refeionce  (4)).  Therefore,  £ 
divided  by  the  energy  density  cannot  be  identically  equal  to  the  surface  wave  group 
velocity;  which  is  a  vector  independent  of  z.  It  will  be  shown,  however,  that  the  group 
velocity'  of  the  surface  wave  is  identical  to  the  velocity  of  energy  transport  of  the  surface 
wave  taken  as  a  whole,  i.  e.  ,  the  gradient  in  the  transverse  wave  number  plane  of  the 
frequency  is  equal  to  the  integral  over  z  of  £  divided  by  tho  corresponding  integral  over 
z  of  the  stored  energy  density.  In  a  manner  analogous  to  that  for  plane  waves  in  aniso¬ 
tropic,  homogeneous  media,  the  relation  between  group  velocity  and  energy  transport 
velocity  for  surface  waves  should  prove  useful  in  formulating  a  ray  interpretation  for  the 
surface  wave  fields  excited  by  a  source.  ^ 

The  proof  of  the  relation  between  the  group  velocity  and  the  eneigy  transport 
velocity  is  furnished  for  two  cor  figurations.  In  the  first  section  cf  this  paper,  the  case 
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considered  is  that  of  a  plane- stratified  medium  filling  all  space,  while  the  second  sec^on 
contains  the  proof  for  the  case  of  a  plane- stratified  medium  filling  the  half- space  above  a 
perfectly  conduct  ng  plane  at  z  =  0. 

In  the  third  aection  of  this  paper,  the  relation  between  the  dyadic  surface  imped¬ 
ance  at  z  -  d,  which  represents  a  plane- stratified,  lossless,  anisotropic  medium  filling 
the  region  0  <  z  <  d  and  bounded  at  z  =  0  by  a  perfectly  conducting  plane,  and  the  power 
flow  and  stored  energy  in  this  region  is  considered.  With  the  help  of  the  developments 
of  the  first  section,  it  is  shown  that  the  power  flow  and  stored  energy  in  the  region 
0  <  z  <  d  are  directly  related  to  the  derivatives  of  the  surface  impedance,  with  respect 
to  transverse  wave  numberc  and  frequency,  and  to  the  components  of  the.  r.  f.  magnetic 
field  transverse  to  z  at  z=d.  The  relation  of  power  flow  and  energy  to  the  surface  react¬ 
ance  apply  for  all  frequencies  and  real  transverse  wave  numbers,  not  just  thc*:e  associated 
with  sur'ace  waves.  In  particular,  for  surface  waves  propagating  above  a  dyadic  imped¬ 
ance  plane,  the  power  flow  and  energy  relations  are  shown  to  bo  significant  in  calculating 
the  energy  transport  velocity. 

The  fourth  section  of  this  paper  is  devoted  to  a  discussion  of  the  dyadic  surface 
impedance  representation  of  a  semi-infinite,  plane- stratified,  lossless,  anisotropic 
medium  for  ranges  of  frequency  and  transverse  wave  numbers  for  which  the  fields  in  the 
medium  are  evanescent  at  infinity.  Again  the  power  flow  and  stored  energy  relation 
involving  the  surface  impedance  are  obtained.  The  power  division  between  the  space 
inside  a.td  outside  of  a  surface  wave  guiding  structure  is  determined  in  terms  of  the  dyadic 
surface  impedances  defined  in  the  third  and  fourth  .actions  of  the  paper. 

The  Appendix  treats  briefly  the  dyadic  admittance  representation  of  a  medium 
above  a  perfectly  conducting  nlane.  At  those  values  of  frequency  and  tran  iverse  wave 
numbers  where  the  impedance  formalism  breaks  down,  the  admittance  formalism  may, 
in  general,  still  be  used.  Power  flow  and  stored  energy  relations  in  terms  of  the  surface 
admittance  are  given. 

STRATIFIED  MEDIUM  FILLING  ALL  SPACE 

A  lossless,  anisotropic,  dispersive,  plane- sirati tied  medium  ?s  assumed  to  fill 
all  space.  It  is  uniform  in  the  x  and  y  directions  and  its  interaction  with  a  monochro¬ 
matic  electromagnetm  field  can  be  described  in  terms  of  the  constitutive  parameters  of 
the  medium,  the  dielectric  tensor  e  and  '.he  permeability  tensor  y  .  Since  the  medium 
is  lossless,  c  and  u  are  Hermilian,  ''  ;  and  because  of  the  assumed  uniformity  in  x 


and  y,  they  a-e  independent  of  these  coordinates.  The  tensors  e  and  |i  0re  analytic 
fu.-Mcions  of  the  angular  frequency  i;  and  are  assumed  to  be  continuous  functions  of  z 
except  for  a  possibly  denumerable  number  of  finite  jumps.  The  •*  dependence  of  c  and 
U  is  further  assumed  to  be  such  that  the  medium  supports  surface  .v»s  propagating 
transversely  to  z.  Such  surface  waves  are  solutions  of  the  source-free  Maxwell  equa¬ 
tions  and  have  the  form 

!*(£»  ]St*  '•«)  e(z;  Jct ,  x) 

H(r;  k t_-  «,•)  |  Mz;  kt>  u) 

where  je  and  h  tend  tc  ~,ero  as  |z|  approaches  infinity  As  used  throughout  this  paper 

p=  x  x  +  y  v  is  the  position  vector  transverse  to  z  and  k=xk+yk  is  a  real 
_  —  o  Ao'  r  —  t  —  o  x  o  y 

transverse  wave  vector.  The  vector  amplitudes  e  and  h  are  required  to  be  such  that  the 
electric  and  magnetic  energy  densities,  as  well  as  Re(exh'),  are  integrable  on  the 
infinite  interval  -®><  z  <  *».  Because  the  dependence  of  IS  and  Id  on  x  and  y  is  e  ■*  — t  £, 

ou  '  “ 

the  electric  and  magnetic  energy  densities  and  Ex  H  are  independent  of  p  .  Furthermore, 
since  the  field  components  transverse  to  z,  and  h^,  must  be  continuous  in  z  across 
any  jump  in  e  or  (j  ,  they  must  be  continuous  functions  of  z  For  simplicity,  e  and  h 
are  assumed  to  be  Rms  quantities. 

In  the  absence  of  sources,  Maxwell's  equations  in  a  medium  described  by  e  and 

u  are 

VxH  =  j  u  c  •  E 
V  x  E  --  -  j  ijb  u  ■  H 

where  the  harmonic  time  dependence  e-*  ‘  has  been  suppressed.  SubsMtuting  Eand  H 

from  equation  (1)  into  (2)  results  in  six  linear  homogeneous  equations  in  the  six  unknown 
field  components,  four  ordinary  differential  equations  in  the  variable  z  and  two  algebraic 
eauations.  For  any  particular  medium  that  can  support  surface  waves,  these  six  equa¬ 
tions  will  have  solutions  satisfying  the  cavity-type  boundary  conditions  e  =  h  =  0  at  |z|  = 
a.id  possessing  the  integrability  properties  described  above,  only  for  restricted  values 
of  the  parameters  k^  and  *  that  satisfy  some  functional  relation  of  the  form 


•.vher.o  in  general  Dy  is  a  regular  function  of  k{  and  a,.  Relation  (3)  is  the  surface 
wave  dispersion  relauon  and  determines  the  y„ssilMe  surface  waves  that  can  propagate 

transversely  to  2  in  the  particular  plane-stratified  medium. 

Let  and  x  be  such  as  to  satisfy  the  surface  wave  dispersion  relation  (3)  and 
consider  neighboring  values  kt  +  *t  and  .  ♦  d«.  also  satisfying  (3).  The  fields  of 
chat  surface  wave  propagating  with  wave  vector  j<t  +  dkt  at  the  frequency  ui+  dx  are 
given,  to  first  order  in  differential  quantities,  by 

-t  +  d-t*  ,;+  da,)  T  h,>  f)  +  5E(r;  k{,  ,) 

H(r;  kt  +  dkt,  x+  d-J  =  H  (r;  kt,  x)  +  5H(r;  k{,  ,t)  (4) 

where  the  variation  6  symbolizes  the  differential  operation 

6  =  dkt  •  V^t  +  dx  £  ,  (5) 

W!th  7kt=^o  *T  +  20  aT;  ■  and  the  Partial  derivatives  of  E  and  H  are  evaluated  at 

(~t*  v)-  Slnce  -t  and  Ht  are  continuous  functions  of  z  for  both  sets  of  values  (k  ,  T 
and  <JSt  +  di$  f  ;t+  dx).  it  is  seen  from  (4)  that  the  variations  &E  and  6H  must  Il^.o 
be  continuous  functions  cf  z.  The  differential  equations  that  6E  and  6H  satisfy  can  be 
found  by  applying  the  variation  5  to  Maxwell's  equations  Recalling  (5)  and  because  € 
and  u  do  not  depend  on  k{,  the  variation  on  Maxwell's  equations  results  in 

*  c  e-  , 

v  x  5H  =  jdx  '  L  +  jxc  •  5E  j 


7  x  =  jd Ji  -  H  -j.au  •  6H 


Conside-  now  the  identity 
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With  the  help  of  equations  (2)  and  (6),  the  right-hand  side  of  (7)  can  be  rewritten  to  give 
the  relation 

jj;  3  .1  ^  jJ;  Sjjkt 

7*  (E  x  6H  +  5E  xH  )  =  -  jd  x(E  •  — —  .  Z  +  H  •  ■  H)  [b) 

o  t*  O  X 

when  the  assumption  that  €  and  u  are  Hermitian  is  used  to  write  6E  •  e  •  E  "  as 
%  #  ~  #  ~ 

E  •  €  •  6 E  ar.d  6H  •  |i  •  H  as  H  •  u  •  5H  .  The  first  term  on  the  right-hand  side  of 

(8)  is  twice  the  time  average  electric  energy  density  w  while  the  second  is  twice  the 

/ 2  g  Q\  ® 

time  average  magnetic  energy  density  w^.  '  *  r  *  As  pointed  3ut>  we  and  are 
independent  of  p.  In  terms  of  the  tctal  energy  density  w  =  wg+  w^i  equation  (8)  can  then 
be  written 

V  •  (E*x  5H  +  &E  x  H*)  =  -j2wd  ;  .  (9) 

The  divergence  term  on  the  left-hand  sid;  of  (9)  is  now  evaluated  by  expanding 

fiE  and  5H  and  subsequently  applying  the  9  operator.  With  the  help  of  (5)  and  (1),  it 

« 

is  seen  that  _j  k  .  f 

5E  =  (fie  -  j  dk  t  P  e)  e  -t  - 

-jk,-  P  ’  (10) 

6H  =  (5h  -  j  dk  h)  e  - 

—  -  -t  _  - 

The  variations  on  the  transverse  vector  amplitudes,  lie  .  and  5ot  »  are  con'  '.nuous 
functions  of  z  since  5_Et ,  5Ht>  £t  and  h  ^  are, a  fact  that  will  prova  useful  later. 

Since  e  and  h  are  independent  of  n  ,  using  the  above  relations  one  has 

V  •  (E*  x  5H  +  iEx  H*) 

=  7  •  [ex  5h  +  6e  x  h*  -  j  d  kt  -p  (e  "x  h  +  e  x  h")] 

=  ^  zQ  •  (e  x  6h  +  5e  x  h  )  -  j2dkt'  p  ^(zq  -  s)  -j2dkt.  s 

(11) 

if  #  is 

where  s(z)  =  Ro(e  xh  )  is  the  real  part  of  the  complex  Poynting  vector  ExH  =  e  x  h  . 
It  is  easily  seen  that  the  term  zq-  £  is  independent  of  z,  i.  e.  ,  -^(_zq-  _p)  =  0,  since  in 
a  source-free  region  filled  with  a  lossless  medium,  the  divergence  of  the  real  part  of  the 
Poynting  vector  is  zero  and,  for  the  plane- stratified  medium  under  discussion,  £  is 
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independent  of  so  that  ~  (xq  s)  =  s)  -  0.  +  Thus  with  the  aic  ■*  (11),  equa¬ 

tion  (9)  '-ecomes 

1  ^  ^  & 

j  j  Q  •  (e  x  6h  +  6e  x  h  )  +  d  k  •  £  =  w  d  x  (12) 

In  the  derivation  of  (12),  the  essential  assumptions  used  are  that  the  medium  be 
lossless  and  that  it  be  plane- stratified  so  that  waves  of  the  form  given  in  (1)  satisfy  the 
source-free  Maxwell  equations.  The  assumption  that  (k^,  x)  and  (kt  +  dkt>  x*  d  ju) 
satisfy  the  surface  wave  dispersion  relation  (3)  serves  to  restrict  the  changes  dk,  jid 
d:jj  in  k^  and  j;  to  a  surface  in  k^.  -  i  space,  so  that  5£t  and  5hf  will  be  continuous 
functions  of  z  for  all  -«<  z  <  ®  and  tend  to  zero  as  |  z|  -  ®  . 

Since  the  wave  vector  k^f  d  k^  and  the  frequency  l  +  d  ;;  satisfy  the  surface 
wave  dispersion  relation,  to  first  order  d  ,  =  dkt-  r(kt>  where  ,r(k  )  is  the 
solution  of  the  dispersion  relation  (3).  Using  this  expression  for  d  x ,  and  after  rear¬ 
ranging,  equation  (1  2)  becomes 

}J  x  5h  +  ^xh  )  =  dkt-  ('v7k  x-  s)  .  (13) 

The  terni  or  the  left  of  (13)  does  not  vanish  identically  so  that  in  gereral  £  w?  ij> 

and  hence  in  the  surface  wave  case  V  ;  cannot  be  interpreted  as  a  locai  energy  velocity. 

—  t 

In  order  to  eliminate  tbj  term  on  the  left-hand  side  of  (13)  we  now  integrate  this 
relation  and  obtain 


where 


and 


j  ^  J  z Q  •  (c  x  5h  +  5e  x  h  )  dz  =  dk(-  (W?^  t:  -  S  ) 

■  00 

00 

S  =  J  s  d  z 

-00 

00 

W  =  J  w  d  z  . 

-as 


(14) 

(15) 

(16) 


Because  sz=z^-  £  is  independent  of  z,  as  discussed  above,  and  is  zero  at  jz|  =  ®, 
since  e  and  h  are  zero  there,  j*  is  zero  for  all  values  of  z.  Thus  £,  and  hence  S, 


Bee  use  e,  and  h.  ace  continuous  functions  of  z,  z  •  s  is  a  continuous  function  of  z 
-t  —  t  —  o  — 

;  d  tnerefore  ~  (_zq  •  £ )  car.nct  have  a  de'‘a  function  behavior  at  the  jumps  of  e  or  u  . 
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are  purely  transverse  vectors. 
Recognizing  that 


$  £  g  jjj 

z  -  (e  x  6h  +  6e  x  h  )  -  z  •  (e.  x5h.H6e.xh.)  , 
—  o  —  —  —  —  —  o  —  t  — t  —  t  —  t 


(17) 


and  using  the  fact  that  e^,  h^,  6e^  and  5ht  are  continuous  functions  of  z,  one  has 


r»  'v 

rr—  z  •  (e  x  6h+  }e  xh  ]dz  =  z  ■  (e.x  5h.  +  6e .  x  h .  ) 

J  o  z  o  *  —  —  — o  — t  — t  — t  — t ' 


(18) 


which  vanishes  as  a  consequence  of  the  boundary  conditions  on  e  and  h  at  |  z|  =  ». 
Hence, 

dfct  (WVR  x-  S)  =  0  (19) 

and,  since  S  is  a  purely  transverse  vector  and  dk  is  arbitrary,  it  follows  that 


V.  t  =  S/W  .  (20) 

-t 

Although  the  real  Poynting  vector  £  can  vary  in  magnitude  and  direction  with  z, 
the  total  real  Poynting  vector  5  is  independent  of  z  and  represents  the  total  surface 
wave  power  flow  across  a  strip  normal  to  S  ,  infinite  in  z  and  having  unit  width.  The 
term  W  represents  the  total  stored  energy  of  the  surface  wave  fields  in  an  infinite 
cylinder,  parallel  to  z,  whose  x-y  cross  section  has  unit  area.  Equation  (20)  thus 
states  that  the  group  velocity  of  the  surface  wave,  is  equal  to  the  velocity  of 

energy  transport  S/W  of  the  surface  wave  as  a  whol£.  This  statement  for  the  surface 
waves  in  plane- stratified  media  replaces  the  relation  9^*  r  s/w  for  plane  waves  in 
homogeneous  anisotropic  media,  and  should  be  useful  in  the  ray  interpretation  of  the  sur¬ 
face  wave  fields  excited  by  a  source  in  anisotropic  plane- stratified  media. 


STRATIFIED  MEDIUM  ABOVE  A  PERFECTLY  CONDUCTING  PLANE 

The  plane- stratified  medium  described  in  the  first  section  is  now  assumed  to  fill 
the  half- space  above  a  perfectly  conducting  plane  at  z  =0.  Again  we  assume  the  z 
dependence  of  e  and  y  to  be  such  that  surface  waves  of  the  form  given  in  (1)  can 
propagate  transversely  to  the  direction  of  stratification.  The  vector  amplitudes  _e  and  h 
of  these  waves  tend  to  zero  as  z  approaches  infinity  and  satisfy  the  boundary  condition 
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£t  =  0  at  z  =  0.  The  electric  and  magnetic  energy  densities,  as  well  as  Re(^  x  h  ),  are 

now  assumed  to  be  integrable  on  the  semi-infinite  interval  0  <  z  <•.  As  discussed  in 

the  previous  section,  e  t  and  h^  must  be  continuous  functions  of  z.  Solutions  of 

Maxwell's  equations  satisfying  the  above  conditions  occur  only  for  values  of  k^  and  u) 

that  obey  a  surface  wave  dispersion  relation,  D  (k  ,  ju)  =  0,  valid  for  the  semi- infinite 

s  —  t 

medium. 

As  in  the  previous  section,  the  fields  at  two  neighboring  sets  of  values,  (kt>  jj) 

and  (kt  +  dkt  ,  «  +  du),  both  of  which  satisfy  the  dispersion  relation,  are  considered. 

Using  equation  (4),  the  fields  E  and  H  at  { k d  1c ^ ,  tu+dui)  are  found,  to  first  order, 

in  terms  of  the  fields  and  their  derivatives,  with  respect  to  k  ,  k  and  m ,  evaluated  at 

x  y 

(kt.  ;u).  Since  the  variations  ir-  the  fields,  fiE_and  fiH_  ,  in  this  problem  also  satisfy  (6), 
equation  (13)  holds  in  this  case  as  well.  Because  the  term  cn  the  left-hand  side  of  (13) 
is,  in  general,  not  zero,  7^  m  again  cannot  be  interpreted  as  a  local  surface  wave  energy 
velocity.  However,  upon  integration  of  (13)  over  the  interval  0  <  z  <®,  the  left-hand 
side  vanishes  and  m  can  again  be  interpreted  as  the  velocity  of  energy  transport  of 
the  surface  wave  as  aSvhole.  To  see  this,  one  recognizes  that  since  £t“0  at  z"0,  fie^ 
must  also  be  zero  there.  Hence,  ising  equation  (17)  and  recalling  that  e  and  h  are  zero 
at  z=  «,  it  is  seen  that 

co  a 

f  -rp-  z  •  (e  x  5h  +  {e  xh  )dz  =  z  •(e’x6h.!,6e  xh)  =0  .  (21) 

I  oz  —  o  —  —  —  —  —  o  —  t  —  t  —  t  —  t 

o  o 

Defining  « 

W  =  J  wd  z  (22) 

o 

and  » 

S  =  j*  sdz  (23) 

o 

(S  ‘being  a  purely  transverse  vector  since  s  =  0)  the  integration  of  (13)  over  the  interv.i 
0  <  z  <  ®  gives,  in  view  of  (21), 

0  =  dk„-  ( W7  w-  S)  .  (24) 

—  t  “ 

Again,  because  dkt  is  arbitrary,  it  follows  that 


7.  m  =  S/W.  . 
-t 


(25) 
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l  hat  is  to  say.  for  surface  waves  above  a  perfectly  conducting  plane,  the  group  velocity 

i  is  equal  to  the  velocity  of  energy  transport  S/W  of  the  surface  wave  as  a  whole.  1 

-•  t 

SURFACE  IMPEDANCE  AND  POWER  FLOW  RELATIONS 


When  formulating  steady-state  electromagnetic  problems  involving  fields  of  the 
form  given  in  (1)  in  a  lossless,  plane- stratified  anisotropic  medium  above  a  perfectly 
conducting  plane  at  z=0,  it  is  sometimes  profitable  to  represent  the  effect  of  the  structure 
below  a  plane  z  =  d  >  0  on  the  fields  in  the  region  z  >d  by  a  surface  impedance  dyadic 
at  z  -  d.  The  impedance  dyadic  Z  may  then  be  employed  as  an  equivalent  boundary  con¬ 
dition  at  z  =  d  when  solving  for  the  fields  in  the  region  x  >  d.  In  this  section  the  relation 
between  the  derivatives  of  the  impedance  dyadic,  with  respect  to  the  spatial  wave  numbers 
and  k^,  and  the  power  flowing  in  the  region  0  <  z  <  d  will  be  established  and  the 
significance  of  this  relation  for  surface  waves  supported  by  such  an  equivalent  impedance 
plane  will  be  pointed  out.  The  relation  between  dZ/ftoi  and  stored  energy  in  the  region 
9  <  z  <  d  will  also  be  established. 

In  order  to  define  Z  and  to  find  its  relation  to  power  flow  and  stored  energy  in 
the  region  0  <  z  <  d,  consideration  is  first  given  to  the  auxiliary  problem  of  finding  the 
fields  in  this  region  when  !jt  of  the  form  given  in  (1)  is  specified  at  z=d.  Thus,  we  look 
in  the  region  0  <  z  <  d  for  the  solution  of  Maxwell's  equations  that  satisfies  the  boundary 

conditions 

--  0  (26) 

at  z  -  0  and  ■»  ,  •,  » 

JUt  -  k  •  o) 

H.(r.t)  =  hde  -  (27) 

at  z  =  d.  All  values  of  k  and  i,  except  those  at  which  Z  is  singular,  are  considered 
(for  further  discussion,  see  the  Appendix).  No  restrictions  are  placed  on  the  fields  in  the 
region  z  >  d.  In  fact,  the  medium  filling  the  region  above  the  plane  z  =  d  may  be  taken 
to  be  arbitrarily  stratified,  since,  with  _H  rigidly  prescribed  at  z  =  d,  the  medium  does 


*If  a  second  perfectly  conducting  plane  at  z  =  d  >  0  is  present,  it  is  easily  seen  that  (25) 
is  still  valid  for  the  fields  between  the  conducting  planes  if  S  and  W  are  now  taken  as 
d  d 

*  (• 

S  :  fa  d  z  and  W  =  i  wdz.  Thus  for  waves  in  a  parallel  plate  wave  guide  filled  with 
o  o 

a  plane- stratified,  lossless,  anisotropic  medium,  the  group  velocity  is  equal  to  the 
velocity  of  energy  transport. 
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not  affect  the  fields  for  0  <  z  <  d.  The  medium  filling  the  region  0  <  z  <d  is  assumed 
to  be  lossless,  uniform  in  x  an-J  y  and  characterized  by  e  and  u,  which  are  analytic 
functions  of  uu- 

Specification  of  the  above  boundary  conditions  >s  sufficient,  in  general,  to  uniquely 
determine  the  fields,  and  hence  the  power  flow  ard  sto~ed  energy,  in  the  region  0  <  z  <d. 
Solving  this  auxiliary  problem  for  arbitrary  polarizations  of  h^  then  permits  a  unique 
determination  of  the  dyadic  surface  impedance  Z.  Having  determined  Z  from  the 
auxiliary  problem,  one  can  now  solve  for  the  fields  above  the  impedance  plane  z  =  d  in 
t“rms  of  Z  ,  the  excitation  in  the  region  z  >  d,  and  the  Loundary  conditions  at  z  =  «. 

The  requirement  that  H  be  continuous  across  z  =d  now  permits  one  to  uniquely  deter¬ 
mine  the  fields,  and  thus  the  power  flow  and  stored  energy,  in  the  region  0  <  z  <d  in 
terms  of  (H  t‘;  _  ^  +  and  the  given  Z  . 

In  practice,  the  auxiliary  problem  need  be  solved  for  only  two  linearly  independent 
polarizations  cf  h^,  since  the  linearity  of  Maxwell's  equations  permits  the  solution  for 
any  other  polarization  of  to  be  expressed  in  terms  of  those  for  the  two  independent 
polarizations.  Thus,  at  each  of  the  values  of  k  and  to  be  considered,  we  solve  for 
the  fields  in  the  region  0  <  z  t  d  when  h^  takes  ou  two  linearly  independent  polarizations, 
e.  g.  ,  hj  =  xq  and  h^  =XQ •  Having  found  the  fields,  which  wiil  be  of  the  form  given  in 
(1),  for  both  polarizations  of  h^,  Z  may  uniquely  be  defined  by  requiring  that  the  relation 

{*t>z  =  d  -  Z  ^ox»t>z=d  <28> 


be  satisfied  for  both  sets  of  fields.  This  requirement  is  equivalent  to  specifying  four 
inhomogeneous,  linearly  independent  aquations  from  which  the  four  unknown  elements  of 
Z  can  be  found.  If  one  now  wisher  to  solve  for  fields  of  the  form  jiven  in  (1),  in  the 
region  z  >d,  relation  (28)  may  be  used  as  a  boundary  condition  at  z  =  d,  which  will  ensure 
that  the  transverse  fields  connect  continuously  to  valid  fields  in  the  region  0  <  z  <d.  That 
is,  if  e  and  h  in  the  region  z  >  d  are  such  that  (28)  is  satisfied,  then,  taking  h^  as 
(h  )  .,  the  corresponding  e  in  the  region  0  <  z  <  d  will  be  such  that  (e  J 


z  =  d 


z  =  d 

H+- 
z  “  d 

Having  thus  defined  Z ,  we  proceed  to  investigate  the  meaning  of  it 3  derivatives 
with  respect  to  k  ,  k  and  w.  To  this  end,  we  assume  the  field*  in  the  region  0  <  z  <d 

^  y 

are  known  and  h  ,  of  (27)  has  been  selected  such,  that  the  derivatives  of  the  fields  with 
—  d 
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"'7=r' t-  X  ^  “d  "  EqUat‘°n  "2’  be  '"-Ployed  where  dk,  and  d„,  i„ 

~  H  ”bUr7  “*  <1«  I-  valid  for  the  fields  “ 

region  -  x  d  sine,  the  assumption,  used  in  deriving  it  are  also  satisfied  in  the 
present  case  -  see  the  text  after  (12).  The  restriction.  Placed  on  dk  and  d„  in  the 
trs  section  are  not  necessary  in  the  present  discussion  since,  a,  previous, y  mentioned 
solutions  Of  Maxwell's  equations  for  which  e,  and  h^  are  continuous  in  . 

.  U.  A  H  "  d  311  ““  *  <'XC,'Vti"S  th'  5taeul*r  POlo-  o'  2.  ..  discussed 

the  Appendix),  thus  ensuring  that  to,  and  ih  are  continuous  function,  of  a  for  all 

,fi.  •'  Jr  lhe  *”  b°"',d"1  SS  °  aPP"aCh"  1  “1  “  '  approaches  0, 

!*  *  P"m,”‘ble  10  ,We*ra"  I12'  “»  closed  interval  0  <  ,  <  d.  Performing  a,, 

integration  and  using  (17)  yields  the  relation 


Here 


il-V<£tx5M^tx^t\  =  d=  wddlt>-  dV  sd  . 


Wd  -  J 


s ,  =  r 


!  s  d  z  . 


Note  that  since  -^-8  =7.s=0  and  Is  1  -  n  *  -n  t  ,,  „ 

.  .  3z  z  -  and  8 3' 7  =  o  “  V°  for  all  0  <  2  <d  and  hence  S, 

is  a  transverse  vector.  —  d 


Since  E  and  H  have  the  form  given  in  (1),  the  impedance  relation 
rewritten  as 


(28)  may  be 


(-t)z  =  d  '  _  '  (£0X  ^t} 


z  -  d 


Applying  the  variation  6  to  the  above  equation 


gives 


1 S-t,Z  =  d  “  2  •  («o  *  tfl^zrd  +  *2-  <*o  *i!t)«  =  d  ’ 

Usir  -t)z  =  d  from  (32)  and  (6et)2  =  d  from  (33),  it  can  be  verified  that 
•  (^tX  6*t  +  «ltx*t\  =  d  =  -  5Z  '  U0*ht)] 


(34) 


It  is  thus  seen  that  and  can  be  found  knowing  only  Z  and  As 

previously  pointed  out,  if  fields  of  the  form  given  in  (1)  exist  in  the  region  z  >d  and 

satisfy  the  impedance  boundary  condition  at  z  =  d,  there  will  be  unique  fields  in  the 

region  0  <z  <  d  that  satisfy  the  continuity  conditions  (h  )  =  (h  )  ,  and 

z=d'z  =  d 

(e  )  =  (e  )  ,  Because  of  the  continuity  of  h  at  z=d,  the  power  flow  anu 

z  =  d'  z  =  d  1 

stored  energy  associated  with  the  fields  in  the  region  0  <  z  <  d  can  be  calculated  from 

relations  (36)  and  (37)  using  =  Lim  h^,  i.  e.  ,  the  limit  of  ht  as  z  approaches 

z  i  d 

d  from  above. 

Since  the  relations  (36)  and  (37)  hold  for  arbitrary  and  they  are  valid,  in 
particular,  for  values  of  kt  and  that  correspond  to  a  surface  wave.  Thus  relations 
(36)  and  (37),  with  appropriate  values  of  and  ji,  furnish  an  alternative  way  of  cal¬ 
culating  that  portion  of  the  surface  wave  power  flowing  in  the  slab  and  that  portion  of  the 
stored  energy  of  the  surface  wave  which  is  in  the  slal/. 


The  impedance  dyadic  Z  is  anti-Hermitian,  i.  e.  ,  the  matrix  representation  for  Z  has 
the  property  that  the  transpose  conjugate  Z+  is  equal  to  -Z.  This  property  follows 
from  the  facts  that  ez  =  0  for  all  z<  d  and  that  the  fields  are  continuous  as  z  approaches 
d  from  below  so  that  Re(e.  xhf)  .  must  be  zero. 

—f  —  t  Z  =  a 
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The  relation  between  power  flow  and  the  derivatives  of  Z  with  respect  to  k  and 
ky  given  in  (36)  does  not  appear  to  have  been  previously  recognized.  While  the  connection 
between  stored  energy  and  5Z/3(d,  to  the  best  of  our  knowledge,  has  not  been  shown 
explicitly  for  the  case  of  traveling  waves,  the  connection  between  stored  power  and  the 
impedance  matrix  of  a  lossless  junction  is  well  known.  ^  ^ 

The  consistency  of  relations  (36)  and  (3?)  for  surface  waves  with  the  results 
obtained  in  the  second  section  will  now  be  shown.  Consider  a  surface  wave  propagating  in 
a  lossless  plane- stratified  medium  above  a  surface  impedance  plane  at  z  =  d.  The  sur¬ 
face  impedance  Z  is  assumed  to  be  known  and  to  represent  the  effect  of  a  plane- 
stratified,  lossless  medium  above  a  perfectly  conducting  plane  at  z=0.  The  surface 
wave  fields  in  the  region  z  >d  are  assumed  to  be  of  the  form  given  in  (1)  with  k^  and  w 
related  through  the  appropriate  surface  wave  dispersion  relation.  The  surface  wave 
fields  satisfy  (13),  which,  when  integrated  over  the  interval  d<z<®,  yields  the  relation 


’■•l-o  ‘  6-t  +  6-t  X  —  t^z  =  d  =  *^t  ^k  ^  lwd  z  *  J-d  z  ^  •  (38) 

_t  d  d 

Since  e(  and  h{  satisfy  the  impedance  condition  (32),  equation  (34)  holds.  Using  (34) 
and  the  fact  that  for  the  surface  wave  dta*  =  dk  •  9,  uj  ,  the  above  equation  can  be  written 

as  "  dZ 

dV  7k  “  {  Jwdz  -  j-|C^o  x-tJ  •  af  ’  (-2oX^t)]z  =  d  } 


f  I>  ;  # 

=  dk.  •  \  sdz+x  -s-  [( z  xh  )■  •  (z  x  h  )] 

— t  1.  J—  —  o  2  —  o  — t  3k  —  o  —  t  Jz  =  ( 


+  y  T(  z  x  h.  )  •  -rt—  •  ( z  xh  )]  .  }• 

io2  —  o  — t  3k  —  o  — t,Jz  =  d  J 


As  discussed  above,  the  terms  containing  3Z/3k  and  3Z/3k  that  appear  in  (39) 

x  ~  y 

are  equal  to  the  x  and  y  components  of  the  power  flow  below  the  plane  z  =  d.  Further 
more,  the  term  containing  3Z/3u.'  is  equal  to  the  stored  energy,  per  unit  area  in  the  x-y 
plane,  below  the  plane  z  =  d,  namely  W^.  Thus  (39)  may  be  written 


•  V1"  {  Jwdz+  Wd}  =  ^t-  {  J^dz  +^d  } 

t  j  j 
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or,  since  dkt  is  arbitrary, 

CP  00 

\  *  {  Iwdz  +  Wd  }  =  Ji'-dz  +  ^d  ■  (41) 

d  d 

Finally,  from  the  definition  of  and  given  in  (30)  and  (31),  equation  (41)  is  seen  to 
reduce  to 

9  9 

j*wdz  =  Jjidz  »  (42) 

1  o  o 

which  is  precisely  the  relation  found  to  hold  in  the  previous  section  for  surface  waves 
above  a  perfectly  conducting  plane  at  z=0.  Hence  if  a  surface  impedance  boundary  con¬ 
dition  representing  a  plane- stratified  medium  above  a  perfectly  conducting  plane  is  used 
when  solving  for  surface  waves,  the  xesultant  group  velocity  V  ^  is  equal  to  the  energy 

transport  velocity  of  the  entire  surface  wave,  not  to  just  that  portion  of  the  surface  wave 
above  the  impedance  plane. 

SURFACE  IMPEDANCE  FOR  THE  CASE  OF  EVANESCENT  WAVES 

In  the  derivation  of  the  power  flow  and  energy  relations  for  a  surface  impedance 
representing  a  plane™ stratified  medium  above  a  perfectly  conducting  plane,  the  presence 
of  the  conducting  plane  served  to  ensure  that  s^  =  C  and  that  ‘he  stored  energy,  per  unit 
area  in  the  x-y  plane,  and  power  flow  are  finite  in  the  region  0  <z  <d  for  all  possible 
polarizations  of  (h^,^  and  all  real  values  of  kt  and  uj.  Since  the  fields  of  evanescent 
waves  in  a  semi-infinite  plane- stratified  medium  also  possess  these  two  properties,  one 
would  expect  power  flow  and  energy  relations  similar  to  (36)  and  (37)  to  exist  in  this 
case  for  the  surface  impedance  representing  the  semi-infinite  medium. 

Let  a  semi-infinite,  plane- stratified,  lossless,  dispersive,  anisotropic  medium 
fill  the  region  above  the  plane  z  =d.  By  analogy  to  the  case  of  the  medium  above  a  per¬ 
fectly  conducting  plane,  consideration  is  first  given  to  the  auxiliary  problem  of  finding 
those  fields  in  the  region  z  >  d  which  satisfy  the  boundary  condition 

j(iyt  -  kt‘p ) 

*£t\s=d  =  -de  (43) 

at.  z  =  d  and  the  boundary  conditions 
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Lim  (E  ,H)  =  0 


In  addition,  we  require  that  Jsidz  and  j'wdz  exist.  The  term  "evanescent,"  as  used 

d  d 

in  the  rest  of  this  pap?r,  will  refer  to  fields  satisfying  (44)  and  the  foregoing  integral 


requirements.  Evanescent  fields  will  also  have  the  property  that  -  0.  The  auxiliary 
problem  is  to  be  solved  for  all  polarizations  of  so  that  the  surface  impedance  may  be 
defined. 


fn  general,  only  for  limited  regions  in  k^-  m  space  w:’l  the  auxiliary  problem 

have  unique,  non-trivial,  evanescent  solutions  that  satisfy  (43)  for  all  h^  and  thus  permit 

definition  of  the  surface  impedance  Z^.  C  ther  values  of  kt  and  -»•  will  not  be  considered 

for  one  of  two  reasons.  First,  in  media  having  an  appropriate  z  dependence,  non-unique, 

cavity-type,  evnnescert  solutions  satisfying  (43)  with  hd=  0  may  exist  for  points  lying 

on  surfaces  in  k^-  o  cpace.  In  cuch  cases,  Zg  will  have  singularities  on  these  surfaces. 

(Discussion  of  such  point3  and  the  derivation  of  a  surface  admittance  formalism  that  is, 

in  general,  regular  at  cuch  points,  are  analogous  to  those  given  in  the  Appendix  for  a 

medium  of  finite  thickness  above  o  perfectly  conducting  plane.  )  Second,  in  some  regions 

of  k  -  (1)  space,  fields  satisfying  (43)  will  not  be  of  the  evanescent  type,  for  most  or  all 
*  T 

polarizations  of  h^.  Thus,  unless  alternate  boundary  conditions  are  specified  at  z  =  « 
s^ch  as  the  radiation  condition,  the  fields,  and  hence  Z  . cannot  be  uniquely  defined.  Even 
If  boundary  conditions  are  impose  at  z=®  and  if  Zg  is  defined  in  this  case^it  is  no  longer 
anti-Hermitian),  the  associated  fields  do  not  possess  the  integration  properties  necessary 
to  derive  simple  power  flow  and  stored  energy  relations. 

Hence  only  those  regions  in  kj  -  iu  space  in  which  the  auxiliary  problem  has 
unique,  non-trivial  solutions  satisfying  (43)  and  (44)  for  all  h^  k  0  will  be  considered 
here.  The  regions  where  the  auxiliary  problem  can  be  solved,  the  nature  cf  which  depends 
on  the  particular  medium  under  discussion,  are  assumed  to  exist  and  to  form  open  sets, 


i.  e.  ,  not  merely  surfaces,  so  that 
within  these  regions. 


k  ,  k  , 
x’  f 


ana  m 


will  be  continuous,  independent  variables 


^An  example  of  a  region  where  no  evanescent  wt  /es  exist  ?.s  formed  by  the  points  in  and  on 
2  1  2  2 

the  cone  = - (k  +k  )  when  the  medium  being  srvdied  is  free  space.  Outside  this 

eoUo  x  y 

cone  unique,  non-trivial  solutions  of  the  auxiliary  problem  exist  for  all  hd- 
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Thus  restricting  kt  and  u  co  those  regions  where  unique,  non-trivial  solutions 
of  the  auxiliary  problem  are  assumed  to  exist  for  all  polarizations  of  h^,  the  surface 
impedance  dyadic  Z  can  be  defined  for  the  semi-infinite  region.  Since  the  linearity 
of  Maxwell's  equations  permits  the  solutions  for  all  to  be  expressed  as  a  super¬ 
position  of  the  solutions  for  two  linearly  independent  polarizations  of  h^,  we  again  need 

consider  only  two  such  polarizations,  e.  g.  ,  h  ,  =  x  and  h  ,  =y  .  From  the  solutions  of 

—  a  —  o  — a  ■*-  o 

the  auxiliary  problem,  which  will  be  of  the  form  given  ir  (1),  for  these  two  polarizations, 
Z  can  be  found  uniquely  from  the  requirement  that  the  relation 


!e.)  ,  =  Z  ■  (-z  x  h  ) 

—  t  z-d  -8  — o  — t  z=d 


(45) 

f 


be  satisfied  by  the  fields  oi  both  solutions.  ^ 

By  analogy  to  the  liscasaion  given  in  the  previous  section,  Z  may  be  used  as  a 

•>•8 

boundary  condition  at  z  =  d  when  solving  for  the  fields  below  this  plane.  Also,  requiring 
to  he  continuous  at  z  =  d  uniquely  determines  the  fields  above  this  plane  when  the 
fields  below  are  known. 

Assuming  Zg  and  the  fields  in  the  region  z  >d  to  te  known,  equation  (12)  is 
employed  iu  finding  the  power  flow  and  energy  relations  in  this  region.  Equation  (i2)  is 
valid  ior  the  fields  ir.  the  region  z  >d  since  the  assumptions  used  in  deriving  it  are 
satisfied  in  the  present  ca3e  —  see  the  text  after  (12).  The  differential  quantities  dkt  and 
dx  in  the  variation  6  are  arbitrary  and  independent  since  ’  k^  and  w  are  independent 
variables  Integrating  (12)  from  z  =  n  to  z  =  ®  gives 

-j|  («t  x  »2it  +  s  x  i?t  )z  =  d  =  ws  -  d_kt  •  S  o  (46) 

where  * 

Wg  =  J  w  d  z  (47) 

d 


In  the  above  relation,  -z  is  used  instead  of  z  ,  as  was  used  in  (28)  and  (72)  for  the 

—  o  —  o 

n  edium  above  ?  perfectly  conducting  plane,  because  -_u Q  is  the  outward  unit  normal  for 

the  configuration  being  considered.  The  conventicn  of  using  the  outward  unit  normal  in 
defining  the  impedance  is  based  on  the  desire  to  have  the  impedance  matrix  be  positive- 
definite  when  loss  is  present  in  the  structure. 


1? 


and  the  purely  transverse  vector  S  is 

—  s 

cr 

Ss  -  J  sd  1  .  (48) 

d 

In  a  manner  similar  to  that  of  the  previous  section,  the  term  on  the  left-hand  side  of  (46) 
can  be  written  in  terms  ot  6Z  g  if  the  anti-Hermitian  property  of  Z  is  taken  into 
account.  That  Z  is  anti-Hermitian  follows  from  the  fact  that  s  =0  for  evanescent 

—  8  Z 

fields  in  a  lossless,  plane- stratified  medium.  In  terms  of  6Zs>  (46)  becomes 


Since  dk  ,  dk  and  dm  are  all  independent,  the  above  relation  can  hold  only  if 
x  y 


C-o  i  *"-o  X-t^  * 


*  aZ8 

air-  •  '-^o*^ 

X 


and 


*  3Z 

+  V  4  ( -  z  xh.)  •  ■~-s  ■  (~z  x  h  )  ]  ,  =  S 

xo  2  — o  — t  ak  -  o  — t  Jz«d  — ! 


1  *  9Z  „ 

-j  —  [(-  z  x  h  )  ■  ■  (-z  x  h  ) J  ,  r  W 

J  2  — o  — t  3uu  —  e  — t  z  =  d  s 


(50) 

(51) 


The  foregoing  relations  should  be  compared  to  (36)  and  (37).  Note  that  if 

instead  of  -z  had  been  usaa  in  (45),  the  abo,re  relation-  would  contain  an  additional 
—  o 

minus  sign. 

The  concept  cf  a  surface  impedance  to  describe  the  effect  of  the  medium  above 
the  plane  z  =  d  on  the  field;-  below  this  plane  can  be  employed  to  derive  the  dispersion 
relation  for  surface  waves.  The  physical  configuration  to  be  consideicd  here  consists  of 
a  plane- stratified,  lossless,  anisotropic  medium  above  the  plane  z-d>0  and  a  second 
plane-!1' ratified,  lossless,  anisotropic  medium  between  a  perfectly  conducting  plane  at 
z  =  G  and  the  plane  z  -  d.  it  will  be  assumed  that  the  values  of  k^  and  ou  of  interest  are 
such  that  the  medium  above  the  plane  c  =  d  is  representable  in  terms  of  an  anti- 
Hermitian  surface  impedance  Zg  that  satisfies  (41*).  This  restriction  on  k ^  and  w  is 
equivalent  to  the  requirement  that  -he  fields  in  the  region  z  >d  be  of  the  surface  wave 
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type  for  all  (h^..^.  (In  special  cases,  surface  waves  may  exist  when  the  fields  in  the 
region  z  >d  are  of  the  surface  wave  type  for  only  one  polarization  of  (hj  Such 

cases  are  not  included  in  the  present  discussion.  )  The  structure  below  the  plane  z  =  d 
ia  assumed  to  be  represented  by  the  surface  impedance  Z  that  satisfies  (32). 

Since  and  h  t  for  the  surface  waves  are  continuous  functions  of  z,  these 
quantities  must  be  the  same  in  both  (32)  and  (45).  Thus  subtracting  these  two  equations 
gives 

»•  (“) 

which  is  a  homogeneous  set  of  two  linear  equations  in  the  two  unknown  elements  of 

(z  x  h  )  For  non-trivial  solutions  of  (52)  to ‘exist  one  requires  that  det(  Z  +  Z  )  =  0, 

which  gives  the  surface  wave  diepersion  relation  D  (k  ,  uu)  =  0.  At  those  values  of  k. 

3  1  “  l 

and  ;c  which  satisfy  the  surface  wave  dispersion  relation,  (j^x  _d  can  be  found. 

If  the  partial  derivatives  of  Z  and  Zg  with  respect  to  k^,  k^  and  uu  are  calculated,  the 
power  flow  and  stored  energy  can  now  be  found  in  each  region  by  using  (36),  (37),  (50)  and 
(51). 

Thus  it  is  seen  t^at  the  knowledge  of  Z  and  Zg  *or  the  lossless  plane- stratified 
structures  previously  described  is  sufficient  to  find  the  surface  wave  dispersion  relation 
and  the  division  of  power  flow  and  stored  energy  between  the  two  regions.  Also,  this 
procedure  can  be  applied  when  the  structure  below  the  plane  z  -  d  is  a  semi-infinite 
medium  whose  regions  in  k^-uu  space,  where  the  reactive  surface  impe  '  ice  maybe 
defined,  intersect  the  corresponding  regions  for  the  medium  above  the  plane  z  =  d. 
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APPENDIX 


In  order  to  define  tuct  Z  which  represents  a  plane- stratified,  lossless,  dispersive, 
anisotropic  medium  above  a  perfectly  conducting  plane,  the  auxiliary  problem,  with 
boundary  conditions  (26)  and  (27),  was  first  considered.  For  most  values  of  the  param¬ 
eters  iS*  and  t^ie  auxiliary  problem  will  have  ur'que  non- trivial  solutions  of  the  form 
given  in  (1)  for  all  hd^  0-  At  the  remaining  values  of  kt  and  Ui,  which  lie  on  surfaces 
in  k t  -  oi  space,  non-unique  cavity-type  solutions  will  exist  for  h  ^  =  0  and  no  solutions 
satisfying  the  boundary  condition  (27)  will  exist  for  all  h^  ^  0.  Thu  non-unique  solutions 
exist  when  the  plane  z  =  d  corresponds  to  a  magnetic  field  null  and  would  give  infinite 
values  for  some  or  all  of  the  components  of  the  impedc-ice  dyadic  Z .  For  thl6  reason 
suck  values  of  and  cu  are  excluded  from  the  consideration  of  Z  .  Although  the  surface 
impedance  formalism  breaks  down  at  these  values  of  k^  and  cu,  a  surface  admittance 
formalism  will  in  general  remain  valid. 

The  surface  admittance  Y  is  the  inverse  of  Z  ,  when  both  Z  and  its  inverse 
exist,  and  is  regular  at  those  values  of  kt  and  u u  for  which  Z  cannot  be  defined.  In 
studying  the  properties  of  Y,  one  would  consider  the  fields  in  the  region  0  <  z  <d  with 
E^, rather  than  H^,  specified  at  z  =  d.  Thus  to  find  the  admittance,  one  requires  that  Y 
be  such  as  to  satisfy  the  relation 


(-5o 


X±t]z 


=  Y 


(^t)z=d 


(53) 


for  two  field  solutions  in  the  region  0  <  z  <d.  The  two  field  solutions  to  be  used  a^e  those 
satisfying  the  boundary  condition 

j(  cut  -  k  •  o ) 

(Et)z=d  =  «de  (S4) 

at  z  =  d  with  ^  taking  on  two  linearly  independent  forms .  If  (53)  is  satisfied  for  these 
two  field  solutions,  because  of  linearity,  it  will  be  satisfied  by  the  solutions  for  all  pos¬ 
sible  e  Those  values  of  k.  a:;d  cu  for  which  non-trivial  solutions  exist  when  e  ,=  0 
—  d  —  t  —  d 

are  excluded  irom  consideration.  Tr.c  Impedance  formalism  may,  however,  be  used, 
in  genera.1.,  at  such  value -i. 

The  energy  and  power  relations  containing  Y  can  be  derived  from  reasoning 
similar  to  that  used  for  Z  .  They  are 


\ 


and 
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l  A  u  * 

»i<V  r  • 


i  * 

J2CV  ar  • 

.1  * 

‘j2[-t  ‘  IS-  •  ■£tJz  =  , 


(55) 


(56) 
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